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ABSTRACT 

We study the interaction of a supermassive black hole (SMBH) binary and a standard 
radiatively efficient thin accretion disk. We examine steady-state configurations of the 
disk and migrating SMBH system, self-consistently accounting for tidal and viscous 
torques and heating, radiative diffusion limited cooling, gas and radiation pressure, 
and the decay of the binary's orbit. We obtain a "phase diagram" of the system as a 
function of binary parameters, showing regimes in which both the disk structure and 
migration have a different character. Although massive binaries can create a central 
gap in the disk at large radii, the tidal barrier of the secondary causes a significant pile- 
up of gas outside of its orbit, which can lead to the closing of the gap. We find that this 

— 1/2 

spillover occurs at an orbital separation as large as ~ 2OOM7 gravitational radii, 
where M, = lO^MyM© is the total binary mass. If the secondary is less massive than 
^ 10^ M0, then the gap is closed before gravitational waves (GWs) start dominating 
the orbital decay. In this regime, the disk is still strongly perturbed, but the piled-up 
gas continuously overflows as in a porous dam, and crosses inside the secondary's orbit. 
The corresponding migration rate, which we label Type 1.5, is slower than the usual 
limiting cases known as Type I and II migration. Compared to an unperturbed disk, 
the steady-state disk in the overflowing regime is up to several hundred times brighter 
in the optical bands. Surveys such as PanSTARRS or LSST may discover the periodic 
variability of this population of binaries. Our results imply that the circumbinary 
disks around SMBHs can extend to small radii during the last stages of their merger, 
when they are detectable by LISA, and may produce coincident electromagnetic (EM) 
emission similar to active galactic nuclei (AGN). 
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1 INTRODUCTION 
1.1 Overview 

Understanding the coevolution of binaries and tlie accre- 
tion disks in which they are embedded is critical in sev- 
eral fields of astrophysics, including planet formatio n and 
migration (|Goldreich fc Tremaind Il98d: IWardI Il997l') . pat- 
terns in plan etary rings (iGoldreich fc Trem aindllQS^), stel- 
lar binaries (|Shu et al'lll987l : iMcKee fc Os trikcr 2 007|), the 
final parsec phase of black hole binaries (|Begelman et al.l 
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Il980l : lEscala et ani2005l : iLodato et ;alll2009l'). stars and black 
holes (BHs) in active galactic nucle i ( Goodman fc Tanll2004l : 
iMiralda-Escude fc Kollmeieill2005l : lLevinll2007l ). and electro- 
magnet ic (EM) counterparts to gravitational wave (GW) 
event s (|Kocsis et al. I I2OO8I : iHaiman et al]|2009l : ISchnittmanI 
I2OIII ). These systems may also provide observational probes 
to test models of the anomalous visco sity in accretion disks 
(jKocsis et al.ll201ll : lYunes et al.ll201lh . 

Despite the long history of the subject, there are very 
few self-consistent analytic models for the coevolution of bi- 
naries and accretion disks, incorporating the fundamental 
physical effects over the long timescales on which the binary 
separation evolves. The standard Q-mo del of radiatively ef- 
ficient turbulent thin accretion disks (jShakura fc SunvaevI 
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1 19731 ) relates the effective kinematic viscosity of the disk 
to the pressure ly oc ap. The viscous evolution of the disk, 
however, is often modeled wi thout considering the pressur e 
dependence of the viscosity (|Lvnden-Bell fc Pringlelll974h . 
Similarly, models of the gravitational interaction between 
the disk, which describe the launching of spiral density 
waves in the disk that remove angular momentum from 
the binary, also do not account for the tidal heating of 
the disk and the correspon ding feedback on the torques 

l|Goldreich fc TremainjilQSOl '). 

In an accompanying paper (jKocsis et all |2012| . here- 
after Paper I) we derive an analytic steady-state model for 
the coevolution of the disk and the orbital migration of 
the se condary, in which we combine a IShakura fc SunvaevI 
lll973l) disk with t he theory of the binary-disk interaction 
bv lGoldreich fc Trcm ainc (1980l ') self-consistently. In partic- 
ular, we adopt the viscosity prescription for standard thin 
a or /? disks, calculate the sound speed and vertical bal- 
ance including both gas and radiation pressure (pgas and 
Prad), adopt a simple analytic approximation to the angu- 
lar in£nMntmn_excl^^ the binary and the disk 
of lArmitage fc NataraianI (I2OO2I'). conside r the viscous and 
tidal heating of the disk ( Lodato et alj |2009| ). and self- 
consistently account for the feedback on the pressure, viscos- 
ity, scaleheight, and the torque cutoff near the secondary's 
orbit, as well as on the migration rate of the secondary. 
We derive a zimuthally averag e d ana lytic disk models which 



recover the iGoodman fc Tan (|2004l ) solution for arbitrary 
P = Pgaa/(Pgaa +Prad) iu the limit that the secondary mass 
rUs approaches zero. The solution is then generalized for 
larger rUs, i.e. to the case when the disk structure is signif- 
icantly modified by the secondary, over multiple accretion 
timescales. 

In this paper, we explore the implications of the new bi- 
nary-|-disk evolution solutions, found in Paper I, for SMBH 
binary systems. By varying the binary parameters system- 
atically, we explore possible distinct behaviors of the disk- 
secondary system. We recover the two limiting cases known 
previously, and identify a new intermediate phase. Low-mass 
objects perturb the disk only weakly, and the linear den- 
sity perturbations lead to t he extensively studied Type- I 
migration of the secondary (|Goldreich fc Tremaind [l98(l ) . 
For very massive objects, the tidal torque clears a gap in 
the disk, and the viscous radial inflow of the gas pushes the 
object inward (known as Type-II migration). A particular 
subclass of the latter is the so-called secondary-dominated 
Type-II migration, in which the secondary's mass exceeds 
the nearby gas mass, causing the migration to slow down, 
and the surface density outside of the gap to build up, before 
it is a ble to efficiently push the object inward jSver fc Clarke! 
1 19951 ). This assumes that the pile-up is 100% efficient and no 
gas can cross the secondary's orbit. We identify a separate, 
intermediate class of migration. Type- 1.5, in which the gas 
piles up significantly outside of the perturber's orbit, but the 
viscosity increases to the point that in steady-state, the gas 
enters the nonlinear gravitational field of the secondary (i.e. 
its Hill sphere), and is able to flow across its orbit. Not sur- 
prisingly, the corresponding migration rate is signiflcantly 
different from both the Type I and II casesQ 



The non-axisymmetric part of the perturbations may be 



In most previous investigations, the gap opening con- 
ditions were based on the comparison o f viscous and tida l 
torques in a weakly perturbed disk (|Crida et al.l l2006l : 
lArmitagdr2007l : lKocsis et aklbOllh . The orbital decay of the 
binary and the evolution of the disk density proflle is then 
coupled; this coupled evolution has not been followed on the 
long time-scales on which the orbit evolves. Although a cav- 
ity may indeed be opened at large radii, in accretion disks 
where the gas build-up outside the cavity is signiflcant, the 
cavity may close after several accretion timescales. Our so- 
lutions allow us to derive the long-term gap opening and 
closing criteria. 

There are many possible observational implications of 
our flndings for black hole binaries, and perhaps also for 
planetary dynamics. Here we restrict our attention to the 
former context. First, gap-closing makes it more likely that 
GW inspiral events are accompanied by EM emission, since 
the binary is embedded in a gaseous disk, with no cen- 
tral c avity, even at t he la st stages of the merger. Previ- 
ously, iLiu et all (|2003l ') and lMilosavlievic fc PhinnevI (|2005l ') 
have argued that as the GW inspiral accelerates beyond 
the rate at which the gas at the edge of the cavity can 
viscously follow the binary, the binary decouples from the 
disk. Consequently, they argued that luminous AGN (ac- 
tive galactic nuclei) or radio emission from jets are expected 
only after the gas has had time to accrete onto the rem- 
nant. The post-merger delay is between years and decades 
for binaries in the mass range 10^-10^ expected to 
be dete ctable by eX/g^/jVCdl^ (lAmaro-Seoane et al.ll201 j: 
see a lso lTanaka et al. l l2010l : lTanaka fc Menoull2Q10l : IShapirol 
|2010D . However, if the central cavity reflUs before the GW 
emission becomes signiflcant (which we flnd is the case in 
the above mass range, in particular), then the gas can ac- 
crete onto the primary and shine much like a normal bright 
AGN, even during the last stages of the merger, producing 
coincident EM counterparts or precursors to LISA sources. 

Second, gas accumulation outside of t he secondary 
leads to a greatly enhanced surface brightness (|Lodato et al.l 
This may help in searches for EM counterparts to 
more massive (10^-10^ M0) SMBH binaries at lar ger sepa- 
rations, still in the gas-driven stage, e mitting GWs in the 
pulsar timing arra y frequency bands jTanaka et al.l I2OI2I : 
ISesana et al.l 20121 ). The EM spectrum of the disk miss- 
ing the emission at high frequencies if the disk has a gap 
ISver fc Clarkdll995l ). and the bright gap ed ge has a char- 
acter istic ultraviolet-optical- infrared profile ( Lodato et al.l 
I2OO9I ). We investigate the brightening as a function of com- 
ponent masses and separations. If these sources produce cor- 
respondingly bright periodic EM variability on the orbital 
timescale, they can be identified in future time-domain op- 
tical/lR surveys. The statistics of many such sources can 
observationally test the migration and GW inspiral rates 
iHaiman et al.l |2009| ). Combining the predictions for the 
variability timescale, disk brightness, and spectrum offers 



either linear or non-linear depending on the magnitude of 
pile-up. In both cases , we adopt the tidal torque model of 
lArmitage fc NataraianI 1 I2OO2I ). reminiscent of Type-I migration, 
which we use with the azimuthally— averaged self-consistent disk 
profile (sec Paper I). 
^ ,http : //elisa-ngo . org/| 
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new independent tests of the physical models of the accre- 
tion plr^Asira^_df^k-satelli^^ 

iGoodman fc TanI (|2004h and iMcKernan et alJ (|2012l l 
pointed out that supermassive stars or intermediate mass 
black holes (IMBHs) in the range 10^-10^ may form in 
AGN disks. These objects would perturb the accretion disk, 
causing a pile up and gap overflow as described here. Our 
studies imply a slow-down of migration for these objects, 
making these systems longer- lived, and thus increasing the 
likelihood of their detection. 



1.2 Relation to previous works 

Without the aim of completeness, we highlight here the sim- 
ilarities and main differences between our study and some 
related papers in the recent lite rature. 

Based o n the solutions of Zel'dovich fc Raized (|l967l ) 
and lPringlj (ligQlh . llvanov et all (|l999l ) showed that the cir- 
cumbinary disk with a binary evolves in a self-similar way 
on scales much larger than the binary, assuming that the 
initial and outer boundary conditions represent an unper- 
turbed stationary a-disk with a fixed M , and assuming that 
a gap is always present which truncates the disk within the 
secondary and causes the radial gas velocity to b e effec- 
tively zero near the secondary (see also I Rafikovll 20 1 3 l . Here, 
we examine the opposite limiting case in which the radial 
gas velocity is nonnegligible near the secondary due to gap 
over flow and the d isk is approximately in a steady-state. 

IChang et al.l (|2010l ) examined the coevolution of the 
disk and the secondary, and found that the former may 
exhibit a rapid brightening in the GW driven regime as 
the binary shepherds ga s inward before merger (however, 
see iBaruteau et al.l l2012h . They included almost the same 
physics as this study and solved the time dependent equa- 
tions in 1 dimension numerically for a binary with 10^ and 
mass components. Their initial conditions were that 
of an unperturbed disk and focussed on the final GW driven 
regime. In this case, they found no gas overflow across the 

iLiu fc Shapirol (|2010l ) constructed an analytic steady- 
state model, and showed that the disk may brighten signif- 
icantly when an object is placed in the disk (locally by a 
factor of fO* for a mass ratio q — 0.1). However, they have 
considered a constant H/r (here H is the disk scaleheight 
and r is the radius) and a viscosity proflle corresponding 
to an unperturbed disk, and neglected the changes in these 
qua ntities due to t he sec ondary. 

iLodato et al] (|2009l ) considered a disk model whose lo- 
cal physics is very similar to ours, but solved the time evo- 
lution numerically in 1 dimension for a very different choice 
of initial and boundary conditions. In particular, they fo- 
cused on the specific case where the SMBH binary is very 
meissive, and has a "one-time" disk that is much less mas- 
sive, compact (spreading over at most a factor of 10 in radii) 
and is not replenished by accreting new material from large 
radii. They found gas pileup outside the secondary, lead- 
ing to a brightening of the outer disk, and a modified spec- 
trum. Their numerical results provide a useful independent 
reference to qualitatively verify our steady-state radial disk 
surface density and scaleheight profiles. For the particular 
binaries they have considered, the migration rate was greatly 
reduced, such that the binary is not transported to the 



GW-driven regime within a Hubble time, and they showed 
that this poses an obstacle aga inst solving the final parsec 
problem (jBegelman et al.lll980l l. In our paper, we examme 
the steady-state configuration, under the assumption of a 
constant accretion rate in the disk (set by the Eddington 
limit near the primary). This assumes a constant mass sup- 
ply from larger radii, and is therefore ve ry different from 
the "one-time" disk in lLodato et al.l (|2009l l . As a result, we 
reach essentially the opposite conclusions: we find a stronger 
pileup, which increases radiation pressure in the disk and 
stabilizes it against gravitational fragmentation, and yields 
a much faster migration (well within the Hubble time). An- 
other practical difference in our study is that we consider 
a broad range of binary masses, mass ratios and semima- 
jor axes, including the radiation pressure dominated regime, 
and we map out the distinct phases for disk structure and 
migration. 

Most studies on planetary migration neglect the tidal 
heating effect and radiation pressure, and are therefore in- 
applicable for our purposes. We investigate AGN accretion 
disks where radiation effects are more significant. Inter- 
estingly, we find that the disk becomes strongly radiation 
pressure dominated outside the secondary, even in regions 
far from the primary, which, in the absence of the sec- 
ondary, would be gas pressure dominated. Therefore, this 
requires treating the fiuid as comprised of both gas and ra- 
diation. The use of a single equation of state parameter, 
as in most papers in planet ary dynamics, b e come s invalid 
in this regime. We note that iD'Angelo et al.l (|2003l ) did ac- 
count for tidal heating and temperature variations in a two- 
dimensional si mulation, but neglected the effe c ts of radi- 
ation pressure . Paardekooper fc Mellemal l|2006l . |2008| ) and 
iKlev fc Cridal (|2008l ) presented results from numerical sim- 
ulations with radiation, showing that tidal heating and ra- 
diation pressure have a significant effect on the migration of 
planets. 

Criteria for gap opening and closing have been 
investigated ext en sively for protoplanetary disks 
(iLin fc Papaloizoul Il986l: lArtvmowicz fc Lubowl Il994l : 
IWardlll997l : iGrida et all 120061 " Previous numerical studies 
typically neglected the effects of gas build-up outside 
the gap, and did not consider self-consistently the effects 
of the excess viscous and tidal heating of the gas, and 
neglected both radiation pressure and the migration 
of the secondary. However, two- and three-dimensional 
simulations have shown that even if the gap opening 
conditions are satisfied, gas can periodically fiow in along 
non-axisymmetric steams into the gap and accrete onto 
the primary and the planet, particularly if the ratio of 
the distance t o the gap edge to the Hi ll radius is of 
order unitv ([Artvmowicz fc Lubowl Il996l: I Lubow et al.l 



19991: iLubow fc D'Angek 



Havasaki et al 



2007 



MacFadven fc Milosavlievid l2008l : iGuadra et al.1 120091 ') . In 



this paper, we revisit the standard gap opening/closing 
conditions in circumbinary accretion disks, in cluding the 
effects of radiation pressure jKocsis et al.|[201ll ). as well as 
gas build-up, tidal heating, and migration. For simplicity, 
we neglect a possible non-axisymmetric inflow if the gap is 
larger than the Roche lobe and also neglect accretion onto 
the secondary. We also do not model the magneto-rotational 
instability (MR I), which may i nfluen c e the conditi o ns for 
gap opening (|Winters et all l2003l : iNoble et al.1 I2OI2I : 
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Figure 1. Gas pile-up and overflow in a circumbinary accretion 
disk witli component masses M, and rria, binary separation r'a, 
and accretion rate M. We distinguish five distinct radial zones in 
the disk: an inner and an outer far zone where the effects of the 
secondary are neghgible, an interior and an exterior near zone, 
where the tidal effects are significant, and an extended middle 
zone, where the disk is still significantly perturbed. Our model 
breaks down inside the secondary's Hill radius, denoted by rj]- 



IShi et al.l [ioil : iFarris et al]|2012l ). These efTects are likely 
to be important, and should be investigated in the future 
in two- or three-dimensional simulations. 

1.3 Outline and conventions 

This paper is organized as follows. In §[2] we briefly lay out 
the physical model that we adopt. In §|3l we present numeri- 
cal and analytic solutions to the disk and identify its distinct 
physical phases. In § lU we elaborate on the gap opening and 
closing conditions in steady-state. In §(5] we discuss the mi- 
gration of the secondary, comparing the new Type- 1.5 solu- 
tion in a continuously overflowing disk with the Type-I and 
Type-II cases. In § [6l we discuss observable implications, 
including the lightcurve, spectrum, the abundance of AGN 
with periodic variability, as well as GW obser vations with 
LISA /NGO and pulsar timing arrays (PTAs, iHobbs et al.l 
We summarize our main conclusions in § [7] The in- 
terested reader can find details of analytic derivations in 
Paper I. 

Our basic notation for the disk and secondary param- 
eters are depicted in Figure [T] We denote the primary and 
secondary mass with M, and rris, and the mass ratio with 
q = rris/M,. We use geometrical units G = c = 1 and sup- 
press factors of G/(? and G/d'' in conversions between mass, 
length, and time units. We use f = r / {GM, / c?) to label 
the radius in gravitational radii. A subscript "s" refers to 
quantities describing the secondary, while the subscript "0" 
refers to quantities in the unperturbed disk around a single 
compact object. 



2 THERMO-HYDRODYNAMICAL 

INTERACTION BETWEEN A DISK AND A 
SECONDARY 

We examine the evolution of the secondary and an az- 
imuthally and vertically averaged Shakura-Sunyaev disk (i.e. 
axisymmetric one-zone disk) in local thermal equilibrium. 



The equations are laid out in detail, and solved both nu- 
merically and under various approximations, analytically in 
Paper I. Here we only briefly summarize the basic concep- 
tual framework and the features of the solutions; we refer 
the reader to Paper I for the complete set of equations and 
a detailed description. 

2.1 Equations governing the evolution of the disk 
structure and the binary's orbit 

We denote the azimuthally-averaged surface density proflle 
of the disk by E(r), and the radial (bulk) velocity of the gas 
in the disk by Vrir) — — M/[27rrE(r)], which is negative as 
gas accretes toward the primary at r = with an accretion 
rate M . We assume a nearly Keplerian disk and denote the 
angular velocity with SI. 

Assuming that the migration velocity of the secondary, 
Vsr is sufficiently slow, the disk evolves through a sequence of 
steady-state configurations, where M(r) = AI is a constant. 
The angular momentum flow in the disk is driven by the 
viscous torques (T^) and the tidal torque exerted by the 
secondary (Td), while the angular momentum loss of the 
secondary is due to the backreaction of the tidal torque and 
gravitational wave losses 



■ 1 



drTi dr - Tgw , 



(1) 
(2) 



Here 



= ~2TTr'^{drO.)iy'E:^3Tvr^nuE, (3) 
a^Td = 27rrAE , (4) 
32 m? 



-7/2 



5 M. 



(5) 



where we adopt the wid ely-used approximation of 
lArmitage fc NataraianI (|2002l ) for the speciflc tidal torque. 



-i/gV'nVVA'' ifr <r,-rH. 
+ lfq^r^n^rt/A'' if r > + ru ., 



(6) 



where 



max(|r — r'a] , H) . 



(7) 

Equation ((6} breaks down in the region near the secondary's 
orbit, i.e. within its Hill radius or tidal radius, \r — rg] < 
(g/3)^/Va. We excise this region from our calculations 
and instead match the interior and exterior solutions, set- 
ting Tv{rs — th) = Tv{rs + rn). This has an effect similar to 
smoothing the torques insid e the H i ll radius, as done pre - 
viou sly in iLin fc Papaloizoul l| 19861 ). ISver fc Clar"k3 l|l995h . 
and iLodato et all (|2009l ). Here, q = rris/M., iJ < r is 
the scaleheight of the disk, and / is a constant calibrated 
with simulations. We conservatively adopt the low value 
/_2 = f/10~^ — 1 for our numerical solutions, but the 
dependence on /_2 is explicitly computed in the analytic 
solutions. 

Given the viscosity !/(r) and the scale-height H{r), 
Eq. ([T]) can be integrated to find the surface density of 
the disk, E(r), for any location of the secondary, and 
Eq. (I2J gives the inward m igration velocity of the sec- 
ondary (|Liu fc ShapirdbOld '). However i/(r) and H{r) are 
not known apriori, since they depend on S(r) and the 
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mid-plane temperature Tc{r). with the standard Shakura- 
Sunyaev ansatz v(r) — acs{r)H{r) for viscosity or its vari- 
ant where the viscosity is assumed to scale with the gas 
(rather than the total) pressure, v{r) = oiCs{r)H{r)j3{r), 
where /3(r) = p^^^sif) / p{r) (also known as a "/3-disk"). Since 
Q disks are thought to suffer from a thermal instability (but 
see iHirose et ahlboogh . in the solutions discussed in this pa- 
per, we follow earlier work and focus on the latter case (also 
known as a "/3-disk" ) . We derive the temperature profile as- 
suming thermal equilibrium between gas and radiation, ra- 
diative cooling which balances the heating associated with 
viscous dissipation tidal heating. These provide a non-linear 
closed set of equations for the steady-state disk profile. In 
practice, we note that all of the equations are nonlinear lo- 
cal algebraic equations at each radius with the exception 
of the angular momentum flux equation ([1} which is a first- 
order ordinary differential equation for Tvir). We solve these 
equations for specific boundary conditions given below, for 
a fixed value of the secondary orbital radius. Once the disk 
profile has been obtained, Eq. © gives the migration veloc- 
ity of the secondary. 

2.2 Boundary conditions 

We look for steady-state disk solutions with a constant M, 
assuming that the radial profile of the disk relaxes on a 
time-scale shorter than the migration time-scale of the sec- 
ondary. We find the corresponding equilibrium steady-state 
configuration of the disk for each orbital radius of the sec- 
ondary, and then assume that the disk proceeds through 
a sequence of such steady-state configurations as the sec- 
ondary migrates inwards. 

We distinguish two types of inner boundary conditions, 
corresponding to whether or not a gap is assumed to be 
present in the steady-state configuration. 

I. Boundary condition without a cavity. We here assume 
that the gas can continuously overfiow, crossing the sec- 
ondary's orbit, and the surface density is finite (S(r) 7^ 0) 
throughout the disk all the way to the primary's inner- 
most stable circular orbit risco ■ We then adopt the zero- 
torque inner boundary condition usually assumed in accre- 
tion disks onto a s i ngle point ma ss (iNovikov fc Thornelll973l : 
iPenna et ai]|2010l : lTanaka|[201ll : IZhu et al.ll2012l ), i.e. we re- 
quire 

r,(nsco) = 0. (8) 

Starting with this boundary condition, we obtain all prop- 
erties of the disk, including the gas velocity profile Vr{r), 
as well as the radial speed of the secondary Vb^. If the 
disk is strongly perturbecjf], the solution is self-consistent if 
|ur('")| ^ K'srI over a wide range of radii around the sec- 
ondary, so that the disk can relax sufficiently rapidly to 
steady-state. In practice, we assume approximate steady- 
state if Ur(Ars) Xvst: (where the constant A > 1 will be 
introduced below). Otherwise, the secondary outpaces the 
nearby gas inflow, and if the disk is strongly perturbed, we 
assume a gap opens. We discuss gap opening in detail in § [l] 



//. Boundary condition for a truncated disk. If the tidal 
torques dominate over the viscous torques near the sec- 
ondary, gas is expelled from the region near the secondary 
and we assume that a circular cavity forms in the disk. In the 
cases we consider, the secondary mass is sufficiently large 
that the gas piles up significantly outside the cavity. The 
tidal toque acting on the inner edge of the disk has a sharp 
cutoff (Eq. We define a characteristic radius in the disk 
outside the gap, where the tidal torque is exerted on the 
dislfl rg — Xrs and assume that the density enhancement 
at this radius tracks the inward migration of the secondary 
in a self-similar way, preserving a constant ratio of radii 
A — rg/rs. This requires that the gas velocity at rg satisfies 

"'■(''g) = — "sr ■ (9) 

rs 

Note that A is not specified by hand ab-initio; it is found self- 
consistently in our solutions (see below and in Paper I). This 
condition can be understood intuitively, since the secondary 
cannot "run away" and leave the outer disk behind (if it did, 
it would cease to be able to torque the disk and would have 
to slow down). Likewise, the gap edge cannot be moving 
closer to the secondary (at least not on on time-scales faster 
than the migration timescale; if it did, then the gap would 
close and the steady-state solution would be inconsistent). 
Although with a moving gap, the disk cannot strictly be in 
steady-state near its boundary, we assume M{r) ~ Af at 
r > rg (see discussion below). 

By construction, only one of the above two boundary 
conditions will lead to a self-consistent solution. We specu- 
late that a real time-dependent binary would evolve through 
the sequence of steady-steady solutions we obtain below - 
switching between the case with and without a gap around 
the transition radii that follows from the above. 

We emphasize that the solutions with a gap are 
somewhat similar to those obtained in previous works 
(|Sver fc Clarke! I1995D . and also that our solutions in this 
regime still suffer from a few possible inconsistencies, as will 
be discussed below. However, the main new result in this pa- 
per is the independent overflowing solution, corresponding 
to the first of the two boundary conditions. As we will argue 
below, the assumptions leading to this regime are relatively 
more robust. The uncertainties about the behavior of the 
disk with a gap could affect only our results for when the gap 
closes (as argued below, we took a conservative approach, 
in the sense that the "overflow" regime may be present for 
a wider range of radii than in our flducial models). 



3 DISK STRUCTURE 
3.1 Numerical solutions 

In Figure [51 we present the most relevant physical param- 
eters in our steady-state disk as a function of radius (we 
neglected Tgw for clarity in this figure) for several choices 
of parameters. For details on how these solutions were ob- 
tained in practice, the reader is again referred to Paper I. 
As the figure shows, the secondary acts as a hydro dam and 



i.e. the viscous torque is greatly reduced (increased) in the near 
zone interior (exterior) to the secondary's orbit relative to the 
solitary disk without a secondary 



In our numerical solutions wc define rg as the radius where the 
tidal torque density drops to 10% of its peak value, although our 
results are insensitive to this precise choice. 
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r \M.] r \M.\ r \M.\ 




r [M,] r [M,] r \M.] 

Figure 2. Steady-state tidally heated disks with the secondary located at = {35, 100, lOOOjA/. (left, center, right panels) for various 
mass ratios q = rris/M, = 10~" with n g {1,2,3,4} (red, green, blue, magenta curves) for Mm = IO^'Mq. The five different rows from 
top to bottom show the local scaleheight, surface density, surface and midplane temperature, surface brightness and tidal heating, and 
the gas to total pressure ratio. The magenta curve with q = 10~* is similar to a solitary disk without a secondary. An inner cavity 
forms in the disk for q = 0.1 and 0.01 for r = 1000 Af, on the right panels. The gas continuously overflows for the other cases shown in 
the figure, but the surface density can still drop significantly near the secondary's Hill sphere (marked by filled circles). (For clarity, the 
curves arc not connected within the Hill sphere, where our model becomes invalid.) 

© 0000 RAS, MNRAS 000, 000-000 



Gas pile-up, gap overflow, and Type 1.5 migration: SMBH binaries 7 



causes gas to pile-up (second row) and heat up (third row) 
relative to the unperturbed solution. Interestingly, the disk 
brightens not only in the near zone of the secondary, where 
the tidal effects dominate, but even much farther away (see 
the panels in the 4th row showing the disk brightness rela- 
tive to the tidal heating rate). This is not surprising, since 
viscosity needs to balance the tidal torque at the boundary 
in equilibrium, and pushing the gas past this obstacle re- 
quires a greatly accumulated local gas-mass, even far from 
this region. Our results are in qualitative agreem ent with 
the radial profiles presented in lLodato et al.l ()2009l ). 

Of the 12 cases shown in Figure [2] with — 
(35, 100, 1000)M. and q = (10"\ 10"^ 10"^ lO""*), a gap 
forms only in the two cases with the largest binary separa- 
tions and mass ratios: Vs = lOOOAf, with q = 0.1 and 0.01 
(red and green curves on the right panels). In these cases, 
there is no inner disk interior the secondary's orbit in the 
steady-state configuration (if there was it would accrete onto 
the SMBH without being replenished from the outside). All 
other cases do not have a gaplf] To avoid confusion, we note 
that the profiles are not shown within the Hill sphere of 
the secondary in Fig. [2] Although our working assumptions 
break down, the parameters are expected to transition con- 
tinuously across this region. 

One may identify visually distinct zones in the radial 
profiles of the physical parameters in Fig. (2] At very large 
radii, which we call the far zone, the disk asymptotes to 
the unperturbed state. Interior to a certain radius, in the 
exterior middle zone, the scaleheight starts to deviate ini- 
tially gradually then quite rapidly as the radiation pressure 
becomes significant relative to the gas pressure. Note that 
the transition to a radiation pressure dominated disk (i.e. 
P < 0.5) is around 600A/. without a secondary. However, 
the disk can become radiation pressure dominated outside 
a secondary much farther out as shown by the 5th row (c.f. 
dotted black line showing jS — 0.5). In the near zone of the 
satellite, the tidal effects become significant. For relatively 
massive binaries, the sharp knee in the tidal heating rate, as 
seen in the 4th row, corresponds to the torque cutoff where 
the midplane pressure gradient shifts the tidal effects out 
of resonance {H ~ I?" — ''sDEI After crossing the radius of 
the secondary's orbit, the accretion velocity is more rapid 
than for an unperturbed disk by many orders of magni- 
tude (see Figure [S] below). The disk profile asymptotes to 
the unperturbed disk profile near the inner boundary close 
to the SMBH. Note however, that the tangential velocity 
{vt ~ (r /M,)~^^^c) and sound speed (cs ~ {H/r)vt) are 
much larger than the radial velocity even in this region. 

We will discuss migration rates in detail in § [5] below, 
but let us mention already which of the 12 cases shown in 
Figure [5] correspond to the well-known Type-I and II cases. 
The two cases = 1000 M. with g = 0.1 or 0.01 have 
a gap, and thus imply Type-II migration. The case with 
q = 10"'' and Ta = 35M, corresponds to Type-I migration, 

^ Here by "gap" we mean regions in which S(r) = (i.e. r < ro 
where > Tb + ru falls outside the Hill radius; see above). Re- 
gardless, note that S(r) can be greatly decreased in the overflow- 
ing solutions near the secondary. 

^ Note that the changes seem deceptively abrupt on a logarithmic 
scale. To see this, consider y = (x — 1)^, a smooth function, that 
exhibits a similar feature on a log— log plot near x = 1. 



where the disk structure lies close to the unperturbed case. 
The 7 cases with q ^ 10"'^ without a gap have a signifi- 
cantly perturbed overfiowing disk, which corresponds to the 
distinct new class of Type 1.5 migration. The remaining two 
cases with q = 10"" and with r^ = (100, 1000)M. are be- 
tween Type I and 1.5. However, comparing the disk-driven 
migration speed with the GW inspiral rate, we find that GW 
emission is significant for many of these binaries, so that the 
steady-state assumption is violated (see Fig.[3]below for the 
conditions under which this occurs). 

3.2 Analytic solutions 

To generalize the numerical solutions given above for arbi- 
trary parameters, we derived in Paper I approximate an- 
alytic formulas. These solutions exist in different regions 
where either the tidal, the viscous torques, or the angular 
momentum fiux is negligible relative to the other two terms 
in Eq. ([l]). In particular, we identify the two far zones, well 
outside and far inside the secondary's orbit, where the effects 
of the secondary are negligible, a single middle zone outside 
the secondary's orbit, where the disk structure is greatly 
modified but where the tidal torque and heating are negligi- 
ble, and two near zones just inside and outside rs, where the 
tidal effects dominate. We distinguish two possible cases for 
the middle zone, depending on whether the disk has a gap 
or if the disk is overfiowing. In the external near zone, we 
likewise have two possible behaviors, depending on whether 
the tidal torque is unsaturated or saturated (i.e. whether 
the torque cutoff is in play: A = |r — rs| or H in Eq. ((6]), 
respectively). We refer the reader to Paper I for detailed 
derivations of the analytic solution for the disk structure in 
each zone; here we provide only the most important results, 
and discuss their implications for SMBH binaries. 

The physical parameters at radius r in the disk are given 

by 

X{r, r„ p) = C a'L\ rht^ r^- fl% q% r% $(r, r„ p) 

(10) 

where X denotes any of \Ti,Tc,ii ,Vr,F,Ty}; r2 and rs2 de- 
note the radial distance from the primary and the orbital 
radius of the secondary in units 100 M, , respectively; and 
^(r, Ts, p) denotes an extra function of the parameters which 
is different in different zones. The constant parameters C 
and Ci and $ are given explicitly in Paper I (see Table 1 
therein). Here we only discuss the structure of the solution. 

The steady-state viscous torque density in different 
regions constitutes the backbone of the analytic solution. 
Every physical quantity in the disk follows directly from 
by simple arithmetic relations. To verify the analytic model 
with the numerical solutions of § 13.11 it is sufficient to ex- 
amine in the various regions. We do this by generating 
the numerical solution for randomly chosen constant model 
parameters (a, rh, f, M,,q, rg), and compare these to the an- 
alytic solution. We find a good agreement between the two 
solutions to within 20% for a wide range of parameters for 
which is substantially modified exterior to the rg. We do 
not develop analytic approximations in the opposite, weakly 
perturbed case, because in that regime other physical pro- 
cesses which have been neglected may be more significant 
(see Type-I migration in § 15. If) . 

The analytic solution is self-consistent in the strongly 
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perturbed case if matches continuously in the different 
zones outside the Hill sphere. In practice, note that T^, in- 
creases in the external near zone quite rapidly with radius, 
both in the tidal torque-saturated and unsaturated regimes. 
At larger radii, the tidal effects vanish (middle zone) and 
is approximately constant. If in the middle zone without 
a gap is lower than with a gap, then this represents the cor- 
rect overflowing solution; otherwise a gap forms. Since for a 
steady-state disk, the migration rate is simply proportional 
to in the middle zone (see Eq. 30 in Paper I), this defini- 
tion is equivalent to that based on the velocities (Eq. [9l). If 
so, the secondary and the gas propagate inward with similar 
velocities to maintain this configuration, while in the former 
case, the secondary moves slower than the inflow speed of 
the gas. We elaborate on these issues related to gap opening 
and migration in § |4] and [5] below. 



We emphasize that the disk parameters can differ dra- 
matically from the unperturbed (or far zone) values not only 
in the near zone, where the tidal effects dominate, but also 
in the middle zone, where tidal effects are already negligi- 
ble. This is analogous to how the water level is raised in 
a dam not just in the immediate vicinity of the dam wall 
but even far away from the boundary, regions in which the 
local hydrodynamics is explicitly independent of the wall. 
The tidal effects of the secondary are also short-range, but 
the corresponding effect is communicated to distant regions 
by setting an effective boundary condition. The middle re- 
gion subtends a large radial range, and is representative of 
the strongly perturbed disk. Here, the viscous torque T™ (in 
all three cases: with a gap, or overflowing with saturated or 
unsaturated torques; see the corresponding solutions T™^, 
T™°=, and T™°" in Table 1 of Paper I) is approximately 
constant with radius as in a decretion disk |Pringlc 199^. 
It is remarkable that the perturbation to the disk can be 
represented by a single number in the middle zone. Here, 
the radial profile of the disk is independent of the details of 
the tidal torque up to this constant factor. 



Rather than using T™, we introduce a more physically 
revealing quantity, k = T^/T^o, representing the fractional 
change of the viscous torque relative to its unperturbed value 
(Ti/O = Mr^rj) in the middle zone. This parameter, which 
we refer to as the dimensionless angular momemntum flux or 
brightening factor, sets the relative brightening of the disk 
and the increase in the scaleheight in the middle zone, as 
well as the migration rate. It can take one of three values, 
depending on whether a gap is open, or whether the gas 
overflows in the saturated torque cutoff regime or in the 



unsaturated regimtQ: 

k^" = 23aJ^mo,l'''M^ 



1/2^-3/8 „_^-3/4 5/8 > -11/16 -7/8 

kr^ = 0.97 al^/^mliM.^/^V^^a'' 9-3" ^fs'"' 



X [-W(-a)] 



13/11 



fcmou 10 — ^■— iTir^/.^j' 
s — 1.6 a_im_iAi^ j_ 



1/2 .5/2 5/2 -1/4 



I-h(^ 



1/3 



-115/24 



-15/2 



(11) 

(12) 
(13) 



Here W(— a) is the Lambert W-function deflned for a > 
1/e = 0.368 approximately as 



\W{-a)\ ^ - In(-a) -f- ln(- ln(a)) 



where 



a = 0.465 a- 1/_ 



-5/4 -13/12 



Mr 



-1/4 5/8 



17/4 



(14) 



(15) 



In practice, 1 < |W(— a)| < 10 for a wide range of parame- 
ters. In the above equations, Sri is the radial distance from r^ 
at which the torque model breaks down near the secondary, 
for which we assume Sri ~ rn. Further, Xr^ is the character- 
istic radius in the near zone outside the gap where most of 
the tidal torque is exerted. In practice, A is of order unity 
given by Eq. (|4ip below. 

The smallest of the three, k^ = mm{k^<^ , k^""" , kf°'^) , 
sets the state of the disk in the middle zone if greater than 
onelfl Note that fc™^ increases quickly with decreasing binary 
separation, r^, but interior to some radius the gap closes and 
the dimensionless angular momentum flux or brightening 
factor is limited by the near-zone torque cutoff, k'^°^ . The 
brightening factor is largest at the gap-closing boundary for 
comparable mass ratios. The torque cutoff does not impose a 
limitation if q is sufficiently small, but the gap can still close 
by shrinking to within the Hill radius (in this case sets 
the brightening factor). 

Figure [3] shows the phase diagram of the disk for dif- 
ferent rs and q, set by Eqs. (|11H13|) for M. = 10^ lO'', 
and 10® M0. Colors indicate the brightening factor, fcs. The 
figure shows large deviations from a solitary disk in many 
cases. The disk brightens in the middle zone by a factor of 
up to 500 or 100 for a binary with M. = 10^ or lO'^M©, 
respectively, which we discuss in § [6] The brightening is less 
significant for binaries with 10® Mq . Note that GW emission 
is neglected here; this modifies the picture at small radii sig- 
nificantly, especially for M, = 10® M0 (see Fig. |4] below) . 

We note that these solutions assume a steady-state. 
However, in Paper I we have shown that the steady-state 
assumption is strongly violated by the inward migration of 
the secondary in the case of a truncated disk with a central 
cavity if M is fixed at the outer boundary (|Sver fc Clarke! 



The brightening factor has a simple radial dependence in the 
middle zone k oc r~^/^. We extrapolate the brightening factor 
in the middle zone to the location of the secondary, typically a 
conservative estimate of the true brightening of the disk outside 
the secondary's orbit as this neglects the excess brightening due 
to tidal heating in the near zone (see Paper I). 
* We do not consider the cases where either fc™^, fc™""^, or 
is less than one. In this case, the analytic solutions need further 
modifications (see Paper I). 
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Figure 3. The phase diagram of the disk for M, = 10^ Mq (top left) 10^ Mq (top right) and 10^ Mq (bottom) and for different binary 
semimajor axes, orbital periods, and mass ratios. Colors show the relative brightening of the disk outside the secondary relative to the 
same region of a solitary disk. Three zones can be distinguished as marked: "gap" in which the tidal torque can support a gap against 
viscosity even in steady-state, "no gap" in which the gap overflows in steady-state with the tidal torque in the linear regime all the way 
to the Hill sphere or "no gap with torque cutoff" in which the torques are saturated outside the Hill sphere. We do not show comparable 
mass ratios {q > 0.1) where the adopted perturbative tidal torque formula and the axisymmetric steady-state approximations are strongly 
violated. 



1 19951 . see however discussion in 15. 2[) . In contrast, steady- 
state is possible in the overflowing case because there ks 
decreases as the secondary moves inwards. A global steady- 
state is possible if ks < fcmax where 



ks 



= 3.2 |7E, 



(5/19 4/19 . -3/19 



M° 



6/19 5/19 -7/19 



1-3 



's2 



(16) 



A similar requirement for a global steady-state in a radiation 
pressure dominated a-disk is less restrictive. For larger ks, 
an approximate steady-state with a fixed M may still hold 
locally in the inner regions of the middle zone, but in this 
case, the outer parts of the middle zone within r ~ fc^rs 
cannot respond to the inward migration of the secondary 
and the steady-state disk model becomes inaccurate there 
(see further discussion in i; 15. 2|) . 

We finally comment on the stability of the accretion 
disk in the middle zone. Although the surface density in- 
creases as gas accumulates outside the secondary as k^^^, 
the radiation pressure grows even more rapidly as p^ad oc k 
in that region (see Eqs. 40-41 in Paper I). The net effect 
is to increase the Toomre-Q parameter as Q oc Cs/E oc 



Prad/S oc k^^^. Therefore, the enhanced radiation pres- 
sure dominates over the increase in surface density, and 
makes the disk more stable against axisymmetric fragmen- 
t at ion relative to an unperturbed disk without a secondary 
(|Tanaka et al.l [20121 ). Although gravitationally stable, the 
disk may become dynamically unstable against global non- 
axisymmetric perturbations if the disk develops a steep pres- 
sure gr adient in the near zone around a high-mass sec - 
ondary (|Papaloizou fc Pringle|[l985l : iGoldreich et al.lll986l ). 
A proper stability analysis, and exploring its implications, 
is left to future work. 



4 GAP OPENING 

Next we elaborate on the gap opening and closing condi- 
tions in more detail. Gap opening is traditionally examined 
by comparing the tidal and viscous torques in an initially un- 
perturbed steady-state thin disk. However, after a sufficient 
amount of gas has built up outside the gap, the enhanced 
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viscous torque may close the gap. We discuss these points 
in turn. 

Gap opening also affects the migration rate of the sec- 
ondary and has several potentially observable signatures 
which we discuss in § [5] and [G] 



4.1 Initial gap-opening 

The velocity of the gas and the secondary can be computed 
starting from Eqs. ([T} and ^ [see Paper I for details]: 



+ 



A 



M 



2 



(17) 



27rrE ' 

drTa dr + To^ ] . (18) 



where negative values represent an inward motion. 

Consider a secondary placed in an initially unperturbed 
disk, and examine the conditions for the local gas flow to be 
reversed. This is equivalent to finding the region Vriro) ^ 
in Eq. (|17|l . which is where the tidal torques dominate over 
the viscous torques, or where the tidal timescale is smaller 
than the viscous timescale in an unperturbed disk. Substi- 
tuting Ti, from Eq. Q and A from Eq. ((G}, and assuming 
that r = Ts -f A where A ^ rs, setting 1;^ = in Eq. (|17p . 
we find 

1/3 

(19) 



Ao 



1^0 



where vo is the unperturbed viscosity near the secondary, 
and fci/o = dlnTi/o/dlnA = 1/2 is the radial exponent of 
the unperturbed viscous torque near r^. Comparing Eq. p9p 
with gas pressure dominated locally isothermal hydrody- 
namical numerical simulations with no tidal heating, give s 
a prefactor in the range 3Ko/f ~ 40-50 (|Crida et al.ll2006D . 
which implies / = 0.03-0.04. 

Two conditions are requirec0 for the tidal torques to 
initially truncate the dis k and create a gap a f ter inserting 
an object into th e disk (|Sver fc Clarke! Il995l : IWardlll997t 
ICrida et al.ll2006D : 



Ao > Ho and Aq > ru ■ 



(20) 



First, the maximum tidal torque, which corresponds to a dis- 
tance Ao ~ Ho, needs to exceed the viscous torque to keep 
gas from fiowing in. Inside this distance, the tidal torque 
saturates due to the "t orque cutoff" and cannot cou nter- 
act the viscous torque (jColdreich fc Tremainel Il980l ) (see 
also discussion below Eq. [5]). Thus, gap opening requires, 
Ao > Ho- Second, the derivation of the torque formula as- 
sumed that the effect of the secondary is a small perturba- 
tion to the gravitational potential; this assumption breaks 
down interior to the Hill sphere. Thus, consistency also re- 
quires Ao > rn. Gas may be expected to accrete onto the 
secondary or to cross the gap, if this condition is violated. In 



addition to violating the conditions in Eq. (|20p . a ga p may 
also close by 3D overfiow if H > Ts. However, Ward (|Wardl 



^ The gap opening conditions are sometimes written as Ao > 



Ho , where rjj 



Ho is the condition for the 



azimuthal perturbations to b ecome nonlinear llWardI Il997l : 
iKorvcanskv &: PapaloizoulllQQSh . We shall not require rj] > Ho 
here over the two conditions in Eq. I I20II . 



1 19861 . [ 19881 ) has shown that using the 2D midplane torque 
gives a result typica lly within 20% to th e vertically averaged 
torque (see however I Jang-Condell fc Sassclov. 20051 ). so that 
the gap closes by midplane inflow in typical cases, and ad- 
ditional criteria are not necessary. 

These gap opening criteria translate into bounds on the 
mass ratio for fixed Ts (|Lin fc Papaloizoulll986l : ICrida et al.l 
I2OO6I) where 



3k.o fHoV^^ , 
a ^ aso = ^/— ( — , and (21) 



(22) 



kvo uo 



If q < Qso the gap is closed by viscosity as the tidal torque 
saturates due to the pressure gradients, and if g < guo it is 
closed by the near-field gravity of the secondary. Thus, after 
releasing an object in a radiation pressur e dominated a or 
P disk, a gap opens if (|Kocsis et al.|[201lh 



qao 



1/2 . 5/2 



>max < 4.5 x 10 



rl/2 5/2 



' 5.7x10-*"-^^-! 



f-2 rl. 



(23) 



fori. in-5 2/5 . 17/10 ,,-1/10 J.-1/2 -29/20 

5/30 -SJ max 2.0 X 10 a_i m_{ M-j ' /_2 r^j , 

A r\ iri-5 4/5 . 2/5 , ,r-l/5 r-l l/io) /r,^\ 

4.0 X 10 a_'i m_'i f_2 r^/, | (24) 

while for gas pressure dominated accretion disks 

g,.,o > max {1.7 x 10-'a]/trh]/,' M,~'^' fZ^^'r'J^ 

4.0 X W-'^a^J^m^J.^M-'^'fiyj'"] . (25) 

Here, the first and second terms in each parentheses corre- 
spond to Qso fEg. I2ip and guo (Eq. [22]), respectively. Note 
that latter is the same for gas and radiation pressure domi- 
nated regimes for a /3 disk because this condition, Eq. (|22p . 
is independent of H, and because f is insensitive to radiation 
pressure in the /3 model. 

Figure |4] shows the gap opening criteria for /3 disks with 
our standard parameters as a function of mass ratio and 
binary separation for M, = 10^'^'^ Mq. The blue solid and 
the dashed curves represent qso and guo, respectively. For 
smaller mass ratios, the disk is practically unperturbed. 

4.2 Steady-state gap-closing 

In the regime where a gap never opens (i.e. weakly perturbed 
disks, q < Qcrito), our solutions do not change the picture. 
However, in the regime when a gap initially opens, the ac- 
cumulation of gas outside the binary changes E, H, v, and 
kv in the gap-opening conditions, Eqs. p9H20p . Let us now 
examine whether the gap can stay open on longer timescales 
once the disk and binary have reached steady-state. 

As the binary migrates inwards the local pressure and 
viscosity change even in an unperturbed disk, which may 
lead to gap opening for mass ratios of order q ~ 10~^ 
within separations Ts ~ 10'' M, (see blue lines in Figure 
However, the usual gap opening conditions Eqs. (|21H22p 
jLin fc Papaloizoul Il986l : ISver fc Clar"k3 Il995l ) may fail in 
the long run and lead to gap-closing for a much wider range 
of masses and radii, if either the scaleheight, the viscosity, 
or the steepness of the radial profile increases during the 
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Figure 4. Ranges of binary separation and mass ratio q for which: (i) a gap can remain open in the disk; (ii) the gap refills after 
gas pileup (labeled as "overflow"), and (iii) a gap cannot open. Different panels show different primary masses (10^, 10^, and 10® Mq) 
as labeled. The type of migration is shown in each region, including the region where GW losses start to become relevant (dotted line) 
and where they dominate (solid black line). Solid and dashed colored lines show the boundary where the tidal torque becomes saturated 
outside the Hill sphere (viscous filling, A = H) and where the gas enters the Hill sphere without becoming saturated (gravitational filling, 
A = rjj), respectively, based on standard criteria using unperturbed disks. The gap (if present) decouples in the GW inspiral regime as 
shown. The disk is gravitationally unstable at large radii as marked. 



pile-up. Note that if only H increased by a factor x while 
V and ku were kept constant, then the mass ratio for gap 
opening would increase by x^^'^ relative to the first terms in 
Eqs. (|23l425fl . Furthermore, if v or changed by a factor 
y, then the the critical mass ratio for gap opening would 
further increase by an extra factor of and y in the first 
and second terms in Eqs. (|23H25[) . respectively. Thus, the 
combination of these effects shifts the solid and dashed blue 
lines in Figure [4] upwards by a factor x'^^'^y^^'^ and y, re- 
spectively. We determine the actual value of x and y in the 
self-consistent steady-state model to derive the necessary 
conditions for gap opening next. 

In a self-consistent decription, gas accumulation outside 
the gap leads to the following effects. First, the gravitational 
torque of a denser disk drives the migration of the secondary 
faster, increasing \vbv\ (see Eos. [4l and I18|) . For a fixed ac- 



cretion rate, M = 2'Kr\vr\T,, the increase of E leads to the 
decrease of 1 1)^1, analogous to the slowdown of the flow veloc- 
ity in a river upstream a dam. The pressure and kinematic 
viscosity increases due to the pile-up, and leads to a smaller 
gap. The tidal heating in the gas outside the disk increases 
the pressure and scaleheight further which can quench the 
tidal torque (see torque cutoff in Eq. [6]). Ultimately, the 
combination of these effects can have several outcomes. 

(i) If the mass of the secondary is very small, the disk is 
weakly perturbed (no gap is present in particular) and the 
object exhibits Type-I migration as usual (see §[5]). 

(ii) The other extreme limit is Type-II migration, where 
l^srl = |ur(Ars)!/A is reached so that the disk and the sec- 
ondary spiral inward in a self-similar way (see Eq. [9|. Here 
the relative gap size, A = rg/rg, is somewhat decreased after 
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the gas has accumulated outside the gap, but it is still much 
larger than the Hill radius. 

(iii) An intermediate possibility is that the gap size 
shrinks over time to within the secondary's Hill sphere as 
gas accumulates outside the gap, and eventually closes. The 
equilibrium steady-state configuration in this continuously 
overflowing case has an increased density and pressure ex- 
terior to the secondary's orbit (see disk profiles in Fig. [2] for 

(iv) Finally, it is also possible that the disk becomes geo- 
metrically so thick {H > r) and/or luminous (L > I/Edd) to 
drive a wind or a three dimensional inflow across the orbit 
of the secondary. 

A necessary condition for the tidal torques to sustain a cav- 
ity in the disk is to satisfy \vsi\ = |?j,.(Ars)|/A, with the 
gap closing if and when |iisr| falls below this value. Here 
Vai, Vr{\ra), aud A represent the final steady-state equilib- 
rium values. This condition can be expressed directly with 
the viscosity in the disk as in Eqs. (|21H22|) . but with the 
viscosity and scaleheight significantly different from their 
original values. We utilize the analytic solution to the com- 
plete disk model (§ 13. 2p as a function of the parameters 
(a, m, M, /, q, r^) and calculate whether these conditions are 
satisfied. In practice, we calculate the corresponding range 
of q and Vs by requiring that the viscous torque at the outer 
boundary of the near zone be equal to that at the inner 
edge of the middle zone with a gap assuming unsaturated 
and saturated tidal torques, respectively (see Paper I for 
details). This is equivalent to finding q or rg for gap open- 
ing at which the dimensionless angular momentum fiux of 
an overflowing solution becomes first equal to that of the 
model with a gap. At the gap opening/closing transition 
k'^'^ = min{fc™°=, in Eqs. iQIMM. Assuming radiation 

pressure dominated /3 disks and negligible GW losses, we 
find that a truncated disk with an inner cavity forms if 



3 ^4/3 r-4/3. 1/3^-1/3 



>4.7x 10-*Q%-'Af7"'Vr2 



-1 ' s2 



^12/31 ^-4/31^3/31 ,,,f-14/31j^j._^^j-104/155 



r,, > 600 M.a_l''fZ*'-='q\^'M^ 



and (26) 

.55 

(27) 



for unsaturated and saturated tidal torques, respectively 
(see Eqs. I14H15I for W and a). Here the first and second 
conditions are analogous to A > rn (or Eq. I22|) and A > H 
(or Eg. I2ip . respectively, but account for the self-consistent 
steady-state radial variations in H, E, and v when calcu- 
lating the viscous and integrated tidal torques for a /3 disk. 
Note that Eq. ^ is defined only for a > 1/e = 0.368, oth- 
erwise the tidal torque is in the unsaturated state. Since a 
depends on qs and rg, Eq. (|27[l is a nonlinear equation for 
ra for a fixed qs- However, the dependence is logarithmically 
weak, typically the last factor in Eq. (|27|) is between 0.3 and 
0.8, making the saturated gap-closing radius vary between 
(200-400)M7 M, for all q for the fiducial disk parameters. 

As mentioned above, for our overfiowing solutions to be 
physically self-consistent, we must also require H < r and 
L < I/Edd, otherwise the radiation pressure would drive an 
outfiow. In § [HI we show that these two conditions are ap- 
proximately equivalent. Assuming there is a gap, a radiation 
pressure dominated /3 disk becomes thick {H > r) for mass 
ratios 



We find this condition to be less restrictive than \vsi\ = 
|ur|/A for unequal masses with q < 0.1. The disk does not 
become geometrically thick if a gap is opened, for our stan- 
dard disk parameters for q < 0.1. The disk thickness is even 
smaller for radiation pressure dominated a disks implying 
that overfiow occurs before the disk drives a wind in this 
case too. These conclusions can be understood by setting 
Ho/rs ~ 1 and uo ~ aCsH/S'' ~ aH^^aP*' ~ ar^VLaP^ in the 
simple gap opening conditions, Eq. H21H22[) . which implies 
gthick <; {'iki^o/fY^'^a^^'^lf^'^. Since the constant of prop or- 
tionality 3fci,o // is of order 40 — 50 in the nonlinear regime 
(|Lin fc Papaloizou|[l98^ : ICrida et al.1120061 ). there is no sec- 
ondary mass that satisfies this cond ition for a > .025 for 
an a disk with 6 = 0. We note that iRafikovl (|2012l ) arrived 
at the opposite conclusion using the same argument but set- 
ting 3k^o/f ~ 1, claiming that the disk typically becomes 
thick and would drive an outfiow before the gap may close 
due to overflow. Furthermore, note that the tidal dissipa- 
tion effects and the steepening of the angular momentum 
flux profile ki, > k^o help further to make the disk prone to 
overfiow before the disk becomes very thick. 

The red line in Fig. [4] shows the boundary where the 
gap remains open in steady-state and where the disk over- 
fiows. Remarkably, the gas overfiows in steady-state and 
closes the gap for q > 10~* for all binary separations be- 
low 1000 M, , and gaps remain open only for much larger 
secondary masses. 

The gap-closing criteria derived above are indepen- 
dent of the initial conditions of the disk assuming that 
the disk is strongly perturbed and relaxes to steady-state. 
However, in Eq. H16|l we have shown that this condition 
does not always hold near the gap-closing boundary and 
a time-dependent study is necessary for a more accurate 
description of gap closing. Furthermore, GW-emission be- 
comes more and more efficient for decreasing binary sepa- 
rations. Figure U shows regions where GW emission starts 
to become relevant and where it dominates. In the region 
marked "gas-fGW driven", the GW-inspiral is slower than 
the unperturbed viscous accretion speed but faster than the 
steady-state gas overfiow. We did not consider this regime 
in detail, but Paper I indicates that there may still be a 
build-up of gas mass, and possibly overfiow here, albeit at 
a slower rate than if GW emission was neglected. In the re- 
gion marked as "GW inspiral" in Fig. |4l steady-state over- 
fiow is prevented by the faster GW-inspiral. In this region, 
the actual state of the disk depends on the initial condi- 
tion. If a gap is open when the binary crosses the crit- 
ical radius where GW-emission becomes dominant, then 
the gap decouples from the binary and remains empty un- 
til the merger. Figure [4] shows that this occurs if the sec- 
ondary mass is larger than approximately ma ^ 4 x 10^ M0 . 
For these masses, an X-ray afterglow is activated once the 
gas refills the gap on the vi scous timescale, typically many 
years after a GW event ([Milosavlievic fc PhinnevI 120051: 



Tanaka et al.l 
2010l : IShapirol 



2OIOI : iLiu fc Shapiro! I2OIOI : I Tanaka &: Menoul 



2010|). However, for smaller secondary masses. 



gthick > 5.4 X 10 -^a 



:1m: 



6/5 
7 



\7/2 3 
A ra2 ■ 



(28) 



Fig. U shows that the gap overfiows and closes before enter- 
ing the GW inspiral phase. We conclude that a gap may 
not be present, and significant X-ray emission accompa- 
nies the GW-emitting stage for these LISA/NGO sources. 
The EM counterpart may be modulated by the GW-driven 
binary jArmitage fc Nataraiaiil |2002| : IChang et all I2OI0I : 
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Farris et al.ll201ll : iBogdanovic et al.ll201ll : iBode et al.ll2012l : 



Noble ct al."2012'; 'Bar uteau et all |2012| : I Giacomazzo et al.l 
2012 : .Farris ct al. 2012l ). but future work is necessary to ex- 
plore this possibility in more detail with initial conditions 
consistent with the overflowing solution. 



5 MIGRATION 

The migration rate of the secondary, v^v, follows from 
Eq. H18|l . where drTd is given in terms of E and H using 
Eq. (|4]-[6| for which we use the steady-state solution derived 
above. The result can be expressed in terms of the dimen- 
sionless angular momentum flux ks in Eqs. (|11H13|) . as 



= -2k, 



Mr, 64 __3 

TTlg 5 



(29) 



where the first and second terms correspond to disk-driven 
migration and to GW-driven inspiral, respectively. 

Figure [5] shows the gas velocity as a function of radius, 
and the migration rate of the secondary for the numerical 
solutions corresponding to Fig. [2] The disk has a gap for 
q = 0.1 and 0.01 when Ts = 10'^ M. , where the nearby gas 
mass and the secondary move radially at a similar speed. 
In this case, the secondary and the disk move inwards to 
maintain a ratio Vr{\rB)/vsi = A, where Arg is the charac- 
teristic radius where the tidal torque is exerted. The increase 
of the gas velocity near the gap edge is due to the fact that 
the steady-state assumption is strongly violated there (see 
further discussion in 15. 2p . At smaller binary separations, 
the gas can flow across the barrier represented by the sec- 
ondary, the disk is no longer truncated, and the migration 
rate is much slower than the gas velocity. 

The gas-secondary interaction has a different character 
when the mass of the secondary is so small that it does not 
perturb the surface density and scaleheight profile signifi- 
cantly (Type-I), and when it is so large that a gap opens 
(Type-II). We first discuss these two limiting cases below, 
and then turn to the intermediate state with a steady-state 
overfiowing disk with a pile-up (which we label "Type-1.5"). 



2011), turbulence and MHD effects (|Nelson fc Papaloizo 



2004 : iLaughlin et all |2004| ), nonlinearities ( Dong ct al 



2011 ) , resonance overlaps (|Rafikov fc Petrovich. 20l3 ) , 



and gas in horse-shoe orbits (Type-Ill migration) 
(|Paardekooper et al.l |2010|). 3D effects for relatively thin 
disks are less significant ( Tanaka et al.|[2002l '). 

We assume that near Tb, E = Ea(r/rs)^^ and H = 
Hs{r /vs)''" , where the exponents 7s and may be slightly 
different inside and outside the secondary's orbit due to the 
tidal effects, y^^u and 73°//, respectively. For unperturbed 
disks, 7_H « and 7s « 3/2 and 7^ ~ —3/5 for a and /3 disks 
respectively in the radiation pressure dominated case, while 
in the gas pressure dominated case 'Jh ~ 21 /20 and 7e ~ 
—3/5 (see Table 1 in Paper I and Eq. (5) in iHaiman et al.l 
I2OO9I ). Note that the scaleheight gradient can be expressed 
in terms of the midplane temperature gradient Tc oc r''^ as 
7H = 3 — 7s + 47T for /3 ^ 1 and 7^- = | + ^jt for /J ~ 1. 

The tidal torque in Eq. (|18[l can be integrated analyt- 
ically in powers of the small quantity H/va in the four do- 
mains of Eq. respectively, where r > or r < (i.e. 
interior vs. exterior zone) and A = _ff or jr — rs| (i.e. satu- 
rated vs. unsaturated torque). To first beyond leading order. 



V,r,I,iu = +|«I.-3 - i^+^-jH] Vl,-2 , (30) 

Wsr,I,is = fl,_3 — ( 2 -f ^ 7h ) W,-2 , (31) 

lie 

«3r,I,cs = — fl,-3 — i -| iH ) "1,-2 , (32) 

Wsr,I,cu = — iwi,-3 — 7h ) «I,-2 , (33) 



where we have expressed the results in terms of 

"I." = 27r/-^E,rfa ( —) . (34) 
M, \ ''s / 

In particular for unperturbed /3-disks, 



5.1 Type-I migration 

Let us first consider a secondary with a mass so small that 
it makes only a small change in the surface density and 
scaleheight profile Eq and Hq and the azimuthal density 
perturbations are in the linear regime (i.e. rn < H). Here 
we estimate the corresponding Type-I migration rate using 
the adopted simple tidal torque model, Eq. (OH For this 
estimate we extrapolate the torque to within the Hill sphere. 

These estimates are included for completeness here, 
but we emphasize that they are very sensitive to 
the approximations used in the tidal torque and the 
disk physics. The true Type-I migration rate (even 
the sign!) can be significantly different for many rea- 
sons i ncluding: thermal effects JPaardekooper fc Mellemal 



I2OO6I: Paardekooper 



( Paardekooper et al.l 



PapaloizorJ 20081 ) . heat diffusion 
bOllI ). inclination (|Bitsch fc Kiev! 



The tidal torque coefficient in Eq. JBJl is / = 
see iGoIdreich fc Tremainj ildSdi ) and Paper I. 



0.80 in this case, 



«l,-2 = < 



A An -1 -4/5 . -7/5 r t fO o 29/ 

440 cms a_i m_^' f-2q--AMj' r^^ 
if ^ < 1 , 



1.2 X 10*cms-^a:i''^m^''f/_2q-3M;'V*2 
if 



^7/5^4/5 



(35) 

Thus, all four domains contribute to the migration rate in 
a nonnegligible way. The total interior and exterior torques 
are 



t^sr,I,i = Wsr,I,iu + l^sr,I,is = 3^!'-^ — (^4 + 7g' — 27)?/^ Vl-2 , 

(36) 



"1,-3 -(7^-275?) 



(37) 



When combining the torques of the interior and exterior 
zones, the leading order term, proportional to 111,-3, drops 
out, irrespective of the exponents. This cancellation makes 
Type-I migration very sensitive to the local disk physics. 
Denoting average (over the inside and outside near zones) 
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quantities by and 7h, 

VsvJ = t^sr,I,i + VsrJ,c ■ 
( 



14 „ 



-VI, -2 



rVj, 



= -(4 + 27S -47H)t;i,_2 (38) 

if ^ < 1, 6 = , 

if /3 < 1, fe = 1 , (39) 

if /3 ~ 1 . 



The last equality corresponds to the exponents in unper- 
turbed a and /3-disks as given above. 

Equation (|38p can be understood as follows. The mi- 
gration rate is a consequence of the opposing repulsive tidal 
effects of the inner and the outer disks. Since the local disk 
meiss is proportional to r'^E oc r^^^^ , the migration rate 
is directed inward for a constant thickness disk (7^ = 0), 
if and only if the local disk mass increases outward. How- 
ever, the scaleheight sets the scale at which the torque is 
suppressed, and acts to reduce the effect. For gas pres- 
sure dominated disks, these estimates match the magni- 
tude of the Type-I migration rate of more accurat e loca l 
torq ue and hydrodynam i cal mo dels of iTanaka et al.l (|2002l ) 
and IPaardekooper et all (|201(l ) if / ~ 0.3 in Eq. © and 
H34p . A different overall offset, mostly due to co-rotation 
torques and pressure gradient effects, directs the migration 
inward in those models for a gas-pressure dominated disk. 
However, recent numerical simulations of non-isothermal op- 
ticall y thick disks with tidal heating find outward migra- 
tion jPaardekooper fc MeUemal |2006| ). and for MHD tur- 
bulent disks_ths_si£n_of_th^ oscillates stochasti- 
cally (|Nelson fc Papaloizoul 120041 '). This highlights the ex- 
treme sensitivity of Type-I migration to subtleties in the 
disk physics for low-mass secondaries. 



5.2 Type-II migration 

Next let us consider the limit that the mass ratio is suffi- 
ciently large (i.e. more comparable masses), the tidal torque 
clears a gap such that the nearby gas accretion velocity 
matches the migration rate of the secondary. If the mass 
of the secondary is larger than the local disk mass then gas 
builds up outside the gap, reducing the accretion velocity 
and increasing the migration rate until the two match. The 
resulting secondary dominated Type-II migration is much 
slower than the accretion velocity of an unperturbed disk. 
Substituting the steady-state viscous torque, T™^ from 



Paper I in Eq. (|18[) we get 

33/4^1/2^1/8 



War, II 



2^1/4(^^p/fc)l/2^1/S 3/S^_ll/160l/Vl/4 



= -550 — a'/,'m'/,'M',/\zl^'X-''^''rU' 



(40) 



assuming GW losses are negligible compared to the gas 
driven migration. Here 1 ^ A < 3 sets the relative distance 
to the nearby gas outside the gap responsible for the tidal 
torque (see Eq. [9]). In Paper I we show that depending on 
whether or not the torque cutoff operates 



A 



1 I n 1 o -4/37 . 4/37 , 



3/37^9/37 22/37 -33/74 



1 + 0.047 a_y^m 



J-2 1-3 



's2 



1/12 ,,fl/6 .1/3 7/12 1/12 



's2 



(41) 

The small blue dots in Figs. [6] and [7| show the migration 
rate of the secondary Vsr and the red circles show local gas 
velocity Vr/X in the numerical solution of §[3Tl]for different 
q and rg. The dotted black curves marked as Type-II in 
Figs. [S] and [7] show that the analytic formula is in excellent 
agreement with the numerical steady-state solution when a 
gap forms. Figure [7] shows that the migration rate is reduced 
relative to the unperturbed disk gas velocity (red curves at 
g — >■ 0) at the given radius. 

Our result for steady-state Type-II migration, Eq. (|40|) . 
is consistent with the secondary d ominated Type-II migra- 
tion rate of ISver fc Clarke! |l99^) , suppressed by a factor 
■^-11/16 < 2 Note that other than M and this weak A de- 
pendence, the Type-II migration rate is insensitive to the 
details of the disk physics and the disk-satellite interaction 
whenever a gap exists. (In particular to the / coefficient in 
the tidal torque Eq. the pressure gradient, the thermal 
state of the disk, and the torque cutoff.) The Type-II mi- 
gration rate is set by the accretion rate, which is set by the 
boundary condition, which we assume is a fixed fraction of 
the Eddington value of the primary. 

Further modifications are possible if the disk is not in 
steady-s tate. The self-s imilar solution of llvanov et~al] (|l999l ) 
(see also lRafiko-Jll2012l ') leads to a slower migration rate (by 
another factor ~ 2 relative to Eq. I40|) . There, the steady- 
state accretion rate is assumed to hold in the far zone outside 
the gap and equal to the rate A^Edd corresponding to the 
Eddington limit near the primary, but the accretion rate 
is reduced in the middle zone relative to the steady-state 
solution, as a fraction of M is continuously used to gradually 
build-up the gas mass outside the gap, which never fully 
reaches the steady-state level. While these models may be 
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Figure 6. Migration rate of the secondary (small blue bullets), radial velocity of the gas (red circles) from the numerical solution 
neglecting GW emission for {M,,q) = (10^ M0,O.l) (left) and (lO'^ Mq, 0.001) (right panel). Analytic migration rates are also shown as 
labeled (Type-II, 1.5, GW) as a function of orbital radius. The magenta line shows the saturated and unsaturated Type-1.5 migration 
rates, black dotted line is the Type-II rate. A gap is present when the blue bullets and red circles overlap. The green dotted line shows 
the GW inspiral rate. The numerical result tracks the analytic Type-1.5 rate when there is no gap. The actual evolutionary sequence of 
the binary is highlighted with a thick yellow line. 
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Figure 7. Same as Fig. |6] but varying q at fixed orbital radii. 



0.1 



more realistic, however, there is no strong reason to believe 
that the far zone accretion rate is related to the Eddington 
limit of the primary, in a scenario with a gap where the gas 
is not present in the vicinity of the primary. If one were to 
assume that M is set by the inner boundary condition in the 
middle zone to a given fraction of MEdd, (implying that this 
may require a correspondingly somewhat larger AI in the 
far-z one), this would have led to the original, quasi-steady- 
state [SyerTciarki (|l995h migration rate. 



Another possibility was considered in iLodato et al.l 

l|2009l ), in which the total amount of gas mass is limited 
to be less than the secondary mass, so that the steady-state 
with Msdd cannot be reached. They argued that the migra- 
tion is much slower in this case, and the disk cannot deliver 
the secondary efficiently to separations where GW emission 



is sufficient to lead to merger in a Hubble time (final parsec 
problem, see §[5]4|. 

5.3 Type-1.5 migration 

In the previous description, Type-I migration requires az- 
imuthally linear perturbations, which is applicable if ru ^ 
H, implying that ? < 3 (H/r)^. By contrast Type-II migra- 
tion assumes a truncated disk with negligible overflow which 
is valid if the gap does not close as discussed in § |4l Let us 
now consider the intermediate Type-1.5 case, when the sec- 
ondary mass is large enough that its tidal torques cause a 
significant gas build-up exterior to the secondary orbit but 
without completely truncating the disk. 

The transition between Type-I and II migration as a 
function of secondary mass was previously investigated ana- 
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lytica lly bv lHourigan &: Ward! (| 19841 ') and lWard fc Hourigan 
1 19891') (see also IWardI Il99'it). numerically by iBate et al 
( 20031 ) ■ and in both ways jy ICrida fc MorbideUil (|2007l ). as- 



suming a constant sound speed and viscosity, in cases where 
the gas pile-up outside the secondary is small. Here we 
discuss Type- 1.5 migration assuming a quasi-steady-state 
derived in Paper 1 in the limit of a large pile-up, self- 
consistently accounting for variations caused by tidal heat- 
ing, which changes the temperature, scale height, viscosity 
in the disk. 

The migration rate is generally proportional to the di- 
mensionless angular momentum flux in the gas driven regime 
(Eq. I29p . This is sensitive to the nearby gas density and 
the characteristic radius at which the torque is exerted, 
VsT <x k cc EA~^, where A ~ max(rH,-H'). If a wide gap 
forms, A is set such that Vsr matches the gas inflow rate 
Var ~ Vr (X M /T, <x u (Type-ll). However, in the overflow- 
ing case, the migration rate becomes slower than Type-H 
as the pile-up is limited by gap overflow, decreasing E and 
H compared to the case with a gap. Relative to an unper- 
turbed disk, although E is decreased close to the secondary, 
it is strongly increased due to pile-up outside of A. How- 
ever, more importantly, the pile-up also increases the vis- 
cous stress and H, which affects the migration rate greatly 
as Vbt oc a H > ru. Accounting for the change in H is 
essential to correctly estimate the migration rate. The com- 
petition of these effects determines the Type-1.5 migration 
rate in the overflowing regime. 

In the self-consistent steady-state overflowing model the 
migration rate can be obtained by substituting ka from 
Eqs. (fT2HT3)) in Eq. (f29)) . In the torque-saturated state, 
which is most relevant for a relatively massive secondary, 
the migration speed is 



^sr , 1.5s 



cm _ 
= 23 — a 

s 



2/11 
1 



j.5/22 ■ 
f-2 1- 



-6/11 

-3 



M7 



23/22 83/44 
^s2 



(42) 



X \W{-a)\ 



13/11 



In the opposite, torque-unsaturated regime for a relatively 
low-mass secondary, 

-35/6 



«sr,1.5u ~ 31 a_ 

s 



2,,r3/2j,5/2 3/2 3/4 



Sri 

Ta 



(43) 



Here GW emission emission is assumed to be negligible, ap- 
propriate in the zone to the right of the dotted line in the 
Vs-q plane in Fig. |4] (i.e. large rs and/or small q). These 
formulae apply only when the tidal torques generate a suf- 
ficiently strong barrier that the surface density of the exte- 
rior disk is greatly modified but a gap does not open. This 
requires ra and q to fall in the zone marked by "overfiow 
(Type-1.5)" in Fig. H 

The migration rate is not always Type-1.5 for other 
choices of q and Tg as shown in Figs. [6]and[7l However, the 
migration rate is approximated to within 20% by Type-1.5, 
Eqs. (|42H43f) . for q > 10^^ and Va 3> r'lsco whenever a gap 
is not present and GW emission is negligible, provided that 
the disk is in a quasi-steady state and the torque model is 
given by Eq. (|6]). These formulae capture the (g, J"s) param- 
eter dependence remarkably well in this range, where the 
Type-I and II rates are inapplicable. Note, however, that 
since our results, Eqs. (|42H43|) . correspond to the case of 
a large pile-up, this is insufficient to investigate very small 



mass ratios where the transition between the Type-I and 1.5 
regimes takes place. 

Figures |6] and [7] show that Type-1.5 migration is gen- 
erally slower than the other velocity scales in the problem. 
Increasing the secondary mass at a fixed orbital radius gen- 
erally increases the amount of gas mass accumulated outside 
the secondary, which causes the gas velocity to be reduced 
for a fixed accretion rate. However, since there is no gap in 
the disk, gas still fiows in faster than the migration rate of 
the object. Type-1.5 migration is slower than the secondary- 
dominated Type-II migration rate, since the gas pileup in 
these overfiowing solutions just outside the orbit is less pro- 
nounced compared to the case with a gap. Finally, it is also 
slower than the Type-I rate if the tidal torques decrease the 
gas density close to the secondary or if the scaleheight is 
increased. 



5.4 GW driven inspiral 

As the separation decreases, the rate of energy and an- 
gular momentum loss associated with G'W emission in- 
creases rapidly. The GW- inspiral spee d (from the standard 
quadrupole formula for GW emission, iPeterslligeA IS 



64 __3 
^qra . 



(44) 



This is shown by a green dotted line in Figs. [6] and [T] Unlike 
Type-I, Type-1.5, or Type-II migration, the GW inspiral ve- 
locity increases toward smaller separations. The GW inspi- 
ral timescale is larger than the Hubble time at large separa- 
tions, where the evolution is driven by the gaseous disk. Gas 
may be responsible for delivering binaries to the separations 
corresponding to the GW-driven regime (but see discussion 
in 15. 2|) . Gas-driven migration may transition to the GW- 
driven regime either in the presence of a gap for comparable 
mass binaries for M > 10^ M0 , or in the overflowing state 
for smaller mass ratios and/or total masses (see Fig. [4|. 

5.5 Summary of migration rates 

We can now summarize and interpret the trends seen for 
different mass ratios and radii in Figs. [S] and [T] At suffi- 
ciently large radii, a gap opens for comparable mass ratio 
binaries. As gas accumulates outside the gap, the radial gas 
velocity is slowed down to match the migration rate of the 
secondary, according to Type-II migration. For smaller sec- 
ondary masses in this regime, the amount of accumulated 
gas is less and the corresponding gas accretion velocity be- 
comes higher and the migration rate faster. At sufficiently 
small secondary mass and binary separation, the tidal torque 
cannot keep the gas outside the Hill sphere and the gap 
closes. The gas accretion velocity in this regime is still far 
slower than the unperturbed value, and the pressure and 
gas density are significantly enhanced outside the orbit. The 
migration here is well approximated by the Type 1.5 rate 
given by Eqs. (|42H43[I . The secondary tends to decrease the 
surface density and pressure within the Hill sphere relative 
to the unperturbed value. The migration rate approaches 
the Type-I rate from below for very small q. Along the in- 
ward migration of the binary, the amount of local gas mass 
torquing the binary decreases. This implies that the mi- 
gration rate decreases with decreasing radius for all three 
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cases (Type-I, 1.5, and II). GW emission dominates over 
gas driven migration in the final stages before merger F'l 



6 OBSERVATIONAL IMPLICATIONS 
6.1 Bolometric luminosity 

The total luminosity of the disk is comprised of the viscous 
and the tidal heating. The energy loss corresponding to in- 
ward migration appears as an extra source of heat in the 
disli, which leads to the brightening of the disk. If the disk 
is radiatively efficient up to an inner radius rmin, then the 
emitted luminosity corresponds to the energy losa^n 



1 GM.M 1 GM.ms , gwn 

^ I mill ' s 



2 Trnij 



(45) 



(46) 



where f = r / {Gc~^ M,) and v^r is the migration rate due to 
gas and GW losses {vst < for inward motion), and is 
the GW inspiral rate neglecting torques from the gas. The 
first term is due to accretion, the second and third terms 
together are due to tidal heating, labeled Lace and ALhcat 
below, respectively. 

Figure [8] shows the two components of the luminosity 
Lace and ALhcat relative to Mc^ as a function of the bi- 
nary orbital time when neglecting GW emission. Here we 
assume that Mc^ — L-Edd = 1.4 x 10''^ M7 erg/sec where 
Mr = M/1O^M0. The point and line styles are the same 
as in Figs. [6] and [7] in the various regimes of migration. The 
figure shows that the secondary can greatly modify the lu- 
minosity of the disk. The accretion luminosity is greatly de- 
creased if a gap forms, but the excess brightness of the outer 
disk, ALheat, can overcompensate for the decrease of Lace. 

The analytic results for the migration rates in §[5]can be 
used to understand the features shown in Figure[Sl If a gap is 
opened and the secondary exhibits Type-II migration, then 
the luminosity associated with accretion scales with the or- 
bital period P as Lace, 11 oc 1/rs oc P~'^l^ and the secondary 
causes an excess ALhcat, 11 oc q^Vsr,ii/r'l oc q^^/^r~^^^^ oc 
gi3/8p-5/4 ^jjg GW-driven regime before gap decou- 
pling, the accretion luminosity is the same but the excess 
tidal heating rate is negligible. At smaller separations the 
accretion luminosity increases, and reaches Afc^/(2risco) as 
the gap closes. As the binary becomes GW-driven at smaller 
separations the gas build-up is decreased outside of the sec- 
ondary and the excess heating is suppressed. This transi- 
tion happens in the region marked as gas-|-GW driven in 
Figs, m at orbital periods P < (81, 14, 10) days for q = 



A separate class of migration for low-mass planets, which 
we did not consider here is Type-Ill migration associated with 
torques exerted by the gas on hors eshoe orbits around the sec- 
ondary jMasset fc Papaloizon|[2003f) . This is significant if the gas 
mass is considerable relative to the secondary mass. Although the 
secondary mass dominates over the local disk mass for BH bina- 
ries embedded in accretion disks of interest here, future studies 
should investigate the effects of the gas crossing the secondary 
orbit (see the recent study bv lRoedig et al]|2012l for the impor- 
tance of torques due to gas streams in the central cavity around 
a SMBH binary). 

We include factors of G and c in this section. 



(0.1,0.01,0.001), respectively, for M. = 10^ M© in partic- 
ular. Once the gap closes, rmin is associated with the ISCO, 
and the accretion luminosity is constant. The tidal heat- 
ing in the overfiowing state is ALhcat, 1.5 — <f "ti 31,1.5/ t'^ oc 
„7/2^-5/4 ^ ^7/2p-5/6 ^j^g unsaturated Type-1.5 state. 



and q 



16/ll_-5/44 



OC gi6/iip-5/66 jj-^ ^jjg saturated case. 



For weakly perturbed disks where the secondary undergoes 



Type-I migration, ALhcat,i = q Vsr,i/rs = q r 



3„9/10 



3p3/5 



in the radiation pressure dominated regime and q r~ = 
g3p-i2/i5 ^j^g pressure dominated regime. 

Figure [8] shows that in some cases ALhcat ~ LEdd- In 
this case, the disk near the secondary becomes brighter than 
the inner accretion disk, so that the disk appears moderately 
super-Eddington. While this may appear contradictory, we 
demonstr ate that this does n ot violat e the local Eddington 
flux limit (|Abramowicz et al.| [l980.: Tanaka fc Men ou 20ld'). 
As shown in § (2] the disk flux can be obtained directly from 
the temperature and the opacity of the disk. If the disk 
is supported vertically by gas and radiation pressure, this 
can be expressed in terms of the disk scaleheight. Assuming 
that the sound speed satisfies H — Cs/fi, the flux is given by 
F = (c/ KjHfi^ {I — j3) (see Paper I). For a nearly Keplerian 
thin disk, the integrated flux from the two faces is then 



47rGM.c 



(1-/?) 



H dr 



(47) 



Here the first factor outside the integral can be identified 
as the Eddington luminosity Lsdd for a source in hydro- 
static equilibrium, while the integral is a geometrical factor 
which is less than unity for a thin disk truncated at the 
ISCO. For a stationary accretion disk without a secondary, 
H{r) is a constant in the most luminous, radiation pres- 
sure dominated regime (/3 « 0), and the integral reduces 
to L = Afc^/(2fmin). More generally, L = r]Mc^ = eL^dd- 
In Fig. [H] we conservatively chose fmin = 6 corresponding 
to a non-spinning BH, and e — rj consistent with AGN 
observations showing t hat e ^ 10%-25% for bright AGNs 
(|Kollmeier et al.ll2006l : iTrump et al.ll2009l '). The luminosity 
increase due the secondary is significant if ff/r is larger near 
the secondary than near the ISCO and if the disk is radia- 
tion pressure dominated. Indeed, Fig. [2]has shown that H/r 
can increase dramatically for q > 10^'^. Eq. ((47| shows that 
a thick, radiation pressure dominated accretion disk can ex- 
ceed the hydrostatic Eddington limit by a logarithmic factor, 
ln(rmax/rmin)- Thls is due to the velocity shear and vorticity 
in the disk, which are neglected in the hydrostatic E ddington 
limit (|Abramowicz et al ] 1 1 980l : I Abramowic j [iooil ') . 

We conclude that the disk luminosity may be modi- 
fied significantly by the orbiting object. Whenever a gap 
forms, the disk becomes either fainter due to the loss of ac- 
cretion onto the primary or brighter due to the accumulated 
gas mass outside the secondary heated by tidal dissipation 
(with the latter effect dominating near the end of the Type 
II migration phase in many cases). The competition of these 
effects generates a unique lightcurve which is sensitive to 
the type of migration as well as to GW losses. These conclu- 
sions assume axisymmetry; an obvious caveat is that non- 
axisymmetr ic accretion streams could generate significan t 
luminosity jHavasaki et al.ll2008l : iRoedig et al.ll201li . I2OI2I '). 
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Figure 8. The bolometric luminosity corresponding to accretion and tidal heating as a function of the orbital period of the secondary 
for (A/,, q) = (10^ Mq, 0.1) and (lO'^ , 10—^) in the various evolutionary stages. The line styles for ALjieat ^^re the same as in Fig. [B] 
black dotted for Type-II migration, solid and dashed magenta for saturated and unsaturated Type-1.5 migration. The cyan line shows 
the total bolometric luminosity Lace + ALheat- Note that our calculations for ALheat neglects GW emission, it breaks down in the 
GW-inspiral regime. Lace is quenched when a gap opens. 
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Figure 9. Disk spectra in units of Lgdcl 10*"^ and 10*^ erg/s for AI, = 10^ and 10^ M, binaries with periods as marked. Different 
curves correspond to different mass ratios as labeled. The dotted curve {q = 10"'^) represents the asymptotic spectrum of a solitary disk. 
The spectrum on the right panel is truncated for q = 0.1 and 0.01 when a gap forms. The optical enhancement at lower frequencies is 
due to the outer disk. 



6.2 Disk spectrum 

We obtain a simple estimate of the disk spectrum by assum- 
ing multicolor black body radiation for each annulus with 
the given surface temperature profile Ts{r), 

and integrating over radius 

= 1^ B[v, Ts{r)] r dr . (49) 

We ignore gravitational redshift and Doppler boost for sim- 
plicity since f S> 1 (however, see discussion below). 

The outer parts of the accretion disk have a much lower 
temperature than the inner parts. Therefore, the brighten- 



ing caused by a secondary aflects the spectrum more promi- 
nently at shorter wavelengths. The disk spectrum can be 
used to disentangle the effects of a secondary from the prop- 
erties of the inner disk. 

Figure [9] shows the disk spectra in two representative 
cases when the orbital period is P = 3 or 110 days for 
M. — 10^ and lO'^M©, respectively. Different curves show 
different mass ratios. The spectrum is truncated at high fre- 
quencies when a gap forms, as the hottest central regions are 
removed from the disk. Relative to an accretion disk with 
no secondary, the disk is much brighter in the ultraviolet, 
optical, and infrared bands for large secondary masses, both 
in the gap forming and in the overflowing cases. The ex- 
tra energy originates from the orbital energy of the binary 
which heats the disk through the gravitational torques on 
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top of the viscous heating. For other binary masses and sep- 
arations assuming the binary is in the gas driven regime, the 
dimensionless angular momentum flux or brightening factor 
is shown in Fig. |3] For smaher binary separations when the 
GW inspiral overtakes the accretion rate, the optical and in- 
frared enhancement associated with the exterior gas pileup 
goes away. 

Note that our estimates neglect accretion onto the sec- 
ondary which can generate an additional spe ctral peak in 
X-rays (jSesana et al.ll2012l : [Tanaka et al.ll2012l 'l. We also ne- 
glected relativistic corrections to the disk spectrum, and 
other features such as spectral lines and the high-frequency 
tail associated with the corona. Since the orbital veloc- 
ity of the secondary on a circular orbit is Vovh = 3 x 
lO^kms"^ rj ^^'^ , the offset and broadening o f spectral lines 
can be quite significant for SMBH binaries jHaiman et an 
20091: IShen fc Loebll2010l : iDotti et al.ll2012l : iMcKernan et all 



2013). Thus, the actual AGN binary spectra shown in Fig. [9] 



may also exhibit relativistically broadened and/or offset 
spectral lines by 6,900 and 9, 500 km s~^ in the left and 
right panels, respectively. The lightcurve is expected to be 
periodically modulated on the binary orbital period (3 and 
110 days respectively in the two panels; see also discussion 
below), with harmonics appearing at ~ 3 and 0.5 times the 
orbital period for mass ratios q > 0.05 (D'Orazio et al. 2012, 
in preparation). These signatures together can help future 
observational efforts to identify SMBH binaries in the over- 
flowing regime. 



6.3 Transient and periodic variability statistics 

We have calculated the rate at which the binary evolves 
toward merger and have shown that the state of the 
disk changes during this process. From this one can pre- 
dict the residence time fres = r/vsi the binary resides 
in the various states. Assuming that the luminosity of 
the accret ion disk is m odulated periodically on the orbital 
timescale jHavasaki et al.i,2007i : iMacFadven fc Milosavlievid 
l2008l : ICuadra et al.ll2009f i. it is possible to predict the rel- 
ative fraction of binari es that exhibit varia bility as a func- 
tion of orbital period (|Haiman et al.l [2OO9I ). One may ex- 
pect to see EM transients when the binary is surrounded by 
gas, and a different distinct population of dimmer or softer 
sources decoupled from the disk assuming that the GW in- 
spiral dominates the evolution inside of a gap in the disk. 
The predictions can be observationally tested using large 
deep variability surveys, such as Pan STARRS or LSST. 

Previouslv. lHaiman et all (|2009l ') has carried out this ex- 
ercise assuming a GW inspiral at short periods and Type-II 
migration for long periods. These predictions have to be re- 
vised in two respects. First, many of the binaries may be in 
the Type-1.5 regime which has a longer residence time. Sec- 
ond, we predict that tidal torques cannot sustain a gap once 
a sufflcient amount of gas has accumulated, which may cause 
gap-refilling, before the binary can "run away" from the gap 
edge in the GW-driven inspiral regime. Without gap-refilling 
or an inner disk around the primary and secondary, peri- 
odic variability is not expected in the GW-driven regime 
with a decoupled outer disk (see however [ Farris et al, I2OIII : 
iBode et all 1 20121 : 1 Noble et al.ll2012l : iFarris et al.ll20ia ). Fu- 
ture studies should investigate the expected variability in 



circumbinary or overflowing disks, where the secondary is in 
the GW-driven regime. 

Figure [TO] shows the residence time the binary spends at 
various orbital periods for different masses and mass ratios. 
Different line styles show the type of migration. The resi- 
dence time increases with period very quickly {trcs oc t^(.^) 
in the GW driven regime, but not so rapidly in the gas driven 
regime. The dependence is flatter for Type-H, and even more 
so for Type-1.5, where it can even decrease with period in 
the torque-cutoff state (at P ~ 1 day for 10^ M© in the right 
panel). This shows that the efficiency at which gas can de- 
liver objects to the GW driven regime is deteriorated by gas 
overflow. The transitions to GW-driven inspiral occurs at 
longer periods for the Type-1.5 than for Type-II migration 
(see Fig. [S| . The fraction of binaries at orbital periods be- 
tween a few days to years is larger if Type-1.5 migration 
operates, relative to the Type-II rate. As a rule of thumb, 
we find this to be the case if one of the black holes is less 
massive than 10^ M0. Conversely, the residence time follows 
the Type-II rate all the way to the GW inspiral regime if 
both objects are more massive than lO'' M0. 

Gap refilling may be witnessed through the increase 
in the bolome tric luminosity and spectral X-ray hardening 
of th e source (|Milosavlievic fc Phinnevll2005l : iTanaka et al.l 
[201^. Our results show that this may occur already be- 
fore merger, partic ularly for 10^ Mq S MBH masses in the 
LISA/NGO range (|Sesana et al.ll2005l ). generating bright 
electromagnetic sources coincident with GWs. Gap refilling 
may occur more rapidly compared to the viscous time with- 
out the pile-up, due to the enhanced gas mass and stress 
outside the gap and the presence of the secondary, implying 
a larger population of birthing quasars for all-sky optical 
(e.g. LSST) or soft X-ray surveys. However, since this oc- 
curs at larger radii where the viscous time is longer, it is 
unclear whether this yields a larger or smaller population 
of birthing quasars with significant brightening during the 
mission lifetime of all-sky optical (e.g. LSST) or soft X-ray 
surveys. 

The significant excess flux in the gas driven regime may 
increase the prospects for identifying SMBH binary sources. 
Future surveys for these sources can provide observational 
evidence for SMBH mergers in the NGO mass range or 
IMBHs orbiting around SMBHs. 



6.4 Gravitational wave measurements 

Gap refilling and Type-1.5 migration has implications for 
GW measurements in several ways. First, we have shown 
that the gap closes much earlier than previously thought. 
This implies that the LISA/NGO binaries may be embedded 
in gas even for nearly comparable mass ratios. This helps the 
prospects of identifying coincident EM counterparts to LISA 
sources. Second, the effects of gas may be identified directly 
from the GW signal itself. The GWs emitted by SMBH bi- 
naries is in the frequency bands detectable by pulsar timing 
arrays (PTAs) if the orbital period is between a few days to 
a few years, and in the planned LISA range if it is between 
a few seconds to a ~half day. Figure [TOl shows that the cor- 
responding mass range is below 10® Mq for LISA and above 
10'' Mq for PTAs. Here we describe how the accretion disk 
affects the GW signal for LISA and PTAs, and discuss the 
implications of gas overfiow and Type-1.5 migration. 
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Figure 10. The amount of time tiie binary spends at various orbital periods during its evolution towards merger. The three panels show 
different mass ratios, and in each panel the three sets of curves show different total masses as marked. For each mass, different line styles 
correspond to different types of migration, GW inspiral (green), saturated and unsaturated Type-1.5 (magenta dashed and solid), and 
Type-II (black dotted). The evolutionary tracks arc highlighted with a thick yellow line. Big black dots mark gap decoupling; for much 
smaller periods, the tidal torques may not generate EM variability. 



6.4.1 PTAs 

The nHz GW backgroun d, measurable by PTAs, i s gener- 
ated by SMBH binaries (|Sesana et al.1 12008| . 120091 ') . which 
emit a stationary signal with angular frequency 2Q. for cir- 
cular sources and many distinct upper harmonics for ec- 
centric sources. The background is affected by the disk 
in two ways: (i) by changing the relative fraction of bi- 
naries at pa rticular orbital radii a ccording to the resi- 
dence time (|Kocsis &: Sesanal I2OIII ). and (ii) by chang- 
ing the eccentricity distribution of sources if a gap is 
opened ([Artymowic 3 1 19921: lArmitage fc NataraianI l2005l : 
ICuadra et al.ll2009l: iRoedig et al.|[201ll . l2012l ). 

Kocsis fc Sesanal ( 20 111 ) calculated the GW background 



for a population of binaries that undergo Type-II migra- 
tion in steady-state a or /3 disks, or the self-s imilar migra- 
tion a ccording to the quasi-stationary model in llvanov et al.l 
l|l999f ). They found that the unresolved GW background 
is typically not reduced significantly in the PTA frequency 
band if the sources undergo secondary dominated Type-II 
migrationl^ This is due to the fact that the background is 
dominated by comparable mass binaries for which (1) the 
slowdown of Type-II migration in the secondary dominated 
regime is more significant and (2) many of these sources 
transition to a GW dominated evolution at large orbital pe- 
riods outside the PTA range. Regarding the individually re- 
solvable GW sources with PTAs, they are typically the most 
massive q ~ 1, ~ 10^ M0 binaries, which are not affected by 
gas in any case at these frequencies. Figure [10] shows that 
Type-1.5 migration in /3-disks is not relevant for PTA mea- 
surements since it does not affect the binaries with masses 
larger than 10® M0 . Since gaps are expected not to close for 
these masses for /3-disks, these sources are expected to ex- 
hibit GW spectra characteristic of eccentric sources. Further 



studies should examine whether gap overflow may occur for 
PTA sources embedded in a disks, and the implications for 
eccentricity. 



6.4.2 LISA/NGO 

As mentioned above, the frequency range of LISA/NGO cor- 
responds to binary orbital periods shorter than a day, the 
sensitivity is best for P ~ 10 min. In this regime, the bina- 
ries are typically already in the GW inspiral regime. Never- 
theless, the effects of gas are imprinted on the frequency 
and phase evolution of the signal as the orb i tal p eriod 
shrinks during the measurement. In lKocsis et al.l |2oTy) and 
lYunes et al.l (|201ll ). we have shown that the corresponding 
GW phase perturbation is very significant (~ lOOOrad/yr) 
relative to the LISA measurement accuracy (~ Irad/yr) 
if the gaseous torque corresponds to secondary dominated 
Type-II migration (particularly for M, = 10® Mq), and less 
pronounced (~ lOrad/yr) but still significant if it corre- 
sponds to Type-I migration. Here we have not examined the 
gaseous torques in the GW-driven regime. In that case, since 
the viscous infiow is slower than the GW-inspiral, the gas 
does not bank-up outside of the secondary^ Therefore, we 
conclude that gas effects ar e expected to be re duced com- 
pared to the Type-II case in lYunes et al.l (|201lh . but future 
studies should examine whether the corresponding phase 
shift is still measurable with the LISA accuracy. 

The GW spectrum may also be affected by the bi- 
nary eccentricity, reminiscent of gas (|Armitage fc NataraianI 
I2OO2I ). We have shown that the disk overflows for typical 
LISA binaries. Further studies should investigate whether 
the eccentricity is excited in the overflowing state, and 



The background could be reduced only if Type-II migration 
occurs on the viscous timescale with no accumulation of gas. 



Gas may bank u p interior to t he orbit in the GW driv en regime 
llChang et al.ll20ld) . see however ISaruteau et al] 1I2OI2I ) and ^ [Ol 
above. 
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whether this remains significant in the GW inspiral regime 
when the binary enters the LISA frequency range. 



7 CONCLUSIONS 

We have examined the self-consistent steady-state structure 
of accretion disks around SMBH binaries, conditions of gap 
closing, and migration rates of the binaries. Our main con- 
clusions can be summarized as follows. 



(i) Gap closing — The secondary represents a strong tidal 
barrier in the accretion disk, which causes a significant ac- 
cumulation of gas in an extended range outside of the sec- 
ondary. Still, the tidal torque cannot counteract the en- 
hanced viscous torque in steady-state for a large range of 
binary masses, particularly for < 10'' M© , and the disk does 
not form a gap. The steady-state exhibits continuous over- 
flow across the orbit. 

(ii) Phase diagram — The disk is described by separate 
formulas in the two main regions in the binary parameter 
space when GW emission is negligible: A gap is present for 
sufficiently large separations and masses (Fig. |3]) but typi- 
cally refills before GW-driven decoupling. We identify fur- 
ther sub-regions in this "phase space" depending on whether 
the tidal torques are saturated. 

(iii) Type 1.5 migration — The gas-driven migration rate 
is Type-I for weakly perturbed disks, Type-II where the gap 
can remain open, and Type-1.5 in the intermediate strongly 
perturbed inffowing state. We have derived analytic formu- 
las for the Type-1.5 rates and have found these to be slow 
compared to both Type I and II rates. 

(iv) optical and infrared excess — The disk fiux is in- 
creased by a factor of up to 100-500 due to gas pileup out- 
side the orbit of the secondary for binaries with M, — 10^- 
10^ Mq in the Type-1.5 and Type-II migration regimes, en- 
hancing the ultraviolet, optical, and near infrared brightness 
of the disk. The enhancement is more pronounced for lower 
M. . The disk can become moderately super-Eddington at 
these frequencies. The excess brightness is mitigated by gas 
overfiow and further suppressed for binaries for which GW- 
emission is significant. This can be used to indirectly indi- 
cate the presence of GWs. 

(v) Periodically variable AGN — The orbital period is be- 
tween 1 day and 10 years for binaries in the Type-1.5 mi- 
gration regime. Periodic variability surveys of AGNs can 
discover these sources and test the predicted spectral sig- 
natures. The statistics of the relative abundance of many 
such sources with different periods can be used to obser- 
vationally test the migration and GW inspiral rate, and to 
estimate the expected LISA merger rate. 

(vi) EM counterparts for LISA — The masses most af- 
fected by gap overflow and Type-1.5 migration is in the 
LISA range (near ~ 10^ M©). The gap closes even for com- 
parable mass ratios before GW emission dominates, imply- 
ing that norm al AGN-like activi ty may be coincident with 
LISA sources (|Kocsis et al.ll2006l '). 

(vii) PTA sources — The gap can remain open for the 
massive 10^-10^ Mq binaries, implying that the GW back- 
ground is not affected by Type-1.5 migration for PTAs. The 
presence of gas can be inferred either from EM surveys or 
indirectly from the GW spectrum, since the gas increases 



the binary eccentricity if a gap is open, which adds orbital 
harmonics to the spectrum. The eccentricity i s expected to 
be much smaller for systems devoid of gas (jSesanal l2010l : 
IPreto et al.ll201ll : iMadigan fc Levinll201 . 

(vii i) Star forma ti on an d final parsec problem — Previ- 
ously, iLodato et al.l (|2009l ) argued that the migration rate 
is greatly reduced by star formation feedback if the disk is 
unstable to gravitational fragmentation and the disk mass 
is smaller than the secondary mass. However, for the low- 
mass, but continuously replenished disks considered here we 
find that the disk is stabilized against gravitational fragmen- 
tation. This is because viscous heating increases the sound 
speed over a large range of radii outside the secondary's or- 
bit, and this effect is more important than the increase of 
surface density in the same region. 

We note that all of these findings correspond to a disk 
model in which the effective viscosity is proportional to the 
gas pressure in the disk (so-called /3-disks). Future stud- 
ies should investigate alternative models in which the vis- 
cosity is proportional to the total gas-|-radiation pressure. 
We also assumed steady-state models, where the accretion 
rate is constant, and set by the Eddington limit of the 
primary. While we expect these assumptions to be justi- 
fied in the overflowing regime, where the gas inflow rate 
is much higher than the migration rate of the secondary, 
the steady-state Eddington accretion assumption is sus- 
pect for circumbin a ry disks with gaps jlvanov et al.lll999l : 
iLodato et al] l2009l : iRafikovl 120121 ). Future studies should 
address comparable mass bina ries where the disk may be 
significantly non-axisymm etric (jMacFadven fc Milosavlievid 
l2008l : ICuadra et all |2009| ) and where the accretion of the 
secondary is non-neglig ible (|Lubow et al.lll999^ ■ and exam- 
ine whether bina ry eccentricity is excited in the overflow- 
ing s t eady-state (lArtvmowicj|l992: lArmitage fc NataraianI 



I2OO5I : lHavasakill2009l : iRoedig et aLlboill . l2012h 
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